Segregation in the annihilation of two-species reaction-diffusion processes on fractal 
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In the reaction-diffusion process A + B —> on random scale-free (SF) networks with the degree 
exponent 7, the particle density decays with time in a power law with an exponent a when initial 
densities of each species are the same. The exponent a is a > 1 for 2 < 7 < 3 and a = 1 for 
7 > 3. Here, we examine the reaction process on fractal SF networks, finding that a < 1 even for 
2 < 7 < 3. This slowly decaying behavior originates from the segregation effect: Fractal SF networks 
contain local hubs, which are repulsive to each other. Those hubs attract particles and accelerate 
the reaction, and then create domains containing the same species of particles. It follows that the 
reaction takes place at the non-hub boundaries between those domains and thus the particle density 
decays slowly. Since many real SF networks are fractal, the segregation effect has to be taken into 
account in the reaction kinetics among heterogeneous particles. 

PACS numbers: 82.20.-w, 89.75.-k, 05.70.1n 
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DifFusion-limited reaction kinetics has been studied for 
fong time as an inter-disciplinary subject. It can be 
a model of electron-hole recombination in semiconduc- 
tors [l| and annihilation of primordial monopoles in the 
early universe 0,01 j etc. The annihilation process involv- 
ing two species A and B of particles A+B is studied 
here, particularly on fractal SF networks. When the den- 
sities of A and B particles arc initially equal, the density 
of each species pA{t) or psit) decays in a power law, that 
is, pA{t) = PB{t) = p{t) ~ In a mean-field approx- 
imation, the density of particles decays as p{t) ^ t~^, 
which is valid when the reaction takes place in Euclidean 
space with the spatial dimension d > dc = 4. For d < dc, 
the exponent a is reduced to a = d/dc, which is less than 
1. The slow decaying behavior of particle density orig- 
inates from the formation of A-rich or _B-rich domains, 
and the reaction takes place at the boundary of those 
domains in Euclidean space 0, S @1 ■ 

In complex networks, however, the particle density p{t) 
can decay faster than the mean- field behavior p{t) ^ t^^ 
in the long time limit Q. This fast decay is caused by 
the existence of hubs at which particles gather through 
the diffusion process and then reaction takes place fre- 
quently. It is noteworthy that the probability of finding 
a random walker at a node is proportional to the de- 
gree of that node [1, Q. For the uncorrelated SF net- 
works, the particle density p{t) was derived analytically 
for A + A ^ and A + B ^ 



III as 
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for 2 < 7 < 3, 
for 7 = 3, 
for 7 > 3, 



(1) 



where 7 is the exponent of the degree distribution 
Pdik) ~ fc-T of the SF networks. 

In this Letter, we demonstrate that when SF networks 
are fractal [l3,[l3|, the segregation of A-rich or S-rich do- 



mains arises, and the particle density decays slowly with 
the exponent a < 1, different from the formula ([1]). Frac- 
tal SF network is a network satisfying the fractal scaling 
Nb{^b) ^ ^s'^^i where Nb is the number of boxes needed 
to cover the entire network with boxes of size Ib- The 
fractal scaling holds when hubs are located separately 
from others in position 14, [l^. Many SF network ob- 
served in real world are fractals. Note that most artificial 
networks including Barabsasi and Albert (BA) model 
are not fractals [17|. In the fractal networks, local hubs 
attract particles and accelerate the reaction. As a result, 
in early time regime, particle density decreases rapidly 
with a > 1. After that period, domains are created in 
which the same species of particles remain, which are the 
majority induced by fluctuations of initial particle densi- 
ties. Then, the reaction takes place only at the boundary 
between those domains, which are not hubs. Thus the 
particle density decays slowly in the long time limit with 
a < 1 . Such segregation behavior can also occur in mod- 
ular SF networks, even if they are non-fractals. Struc- 
tural feature of the modular network, being composed of 
a large number of links within modules but a small num- 
ber of links between modules, hampers the diffusion of 
particles across modules. 

To study the two-species reaction A-\- B ^ ow frac- 
tal SF networks specifically, we first recall the previous 
studies Q of the reaction kinetics taking place on frac- 
tal structure embedded in Euclidean space. In this case, 
the formula p{t) ^ t-d/dc y^^y replaced with 



p{t)^t-''^/\ 



(2) 



where dg is the spectral dimension of the fractal struc- 
ture, ds is related to random walk dimension dw and 
fractal dimension df as dg = 2df/dw The random 
walk dimension is defined through the anomalous power- 
law relationship between the mean-square displacement 
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if- (t ) ) of a diffusing particle and time t as {f'{t)) t^/'^^ . 
The formula ([2]) has been questioned, however, because 
it does not take into account of structural features in a 
given fractal structure such as the degree of ramification. 
Nevertheless, it appears that numerical results are es- 
sentially in agreement with this prediction ^ for many 



cases 
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In this Letter, we show that in contrast 
to the standard random SF network cases, for the fractal 
SF networks we study here, the particle density decays 
in the form given by ([2|). 

Here we first generate a fractal SF tree structure 
through the multiplicative branching process. At each 
branching step, a node creates its m branches (offsprings) 
with probability ~ 'm~^ (jn > 1). It has to satisfy the 
criticality condition (m) = X]m=o '"^Pm = 1 [Si ■ Then, 
the resulting tree structure is a SF tree with the degree 
exponent 7. Such a random critical branching tree struc- 
ture is a fractal SF network with the fractal dimension 
d/ = (7 - l)/(7 - 2) for 2 < 7 < 3 and d/ = 2 for 
7 > 3. The spectral dimension is dg — 2(^ — l)/(27 — 3) 
for 2 < 7 < 3 and 4 = 4/3 for 7 > 3 [Ifl, llH . 

We measure particle density p{t) as a function of time 
t in the form, 
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p{t) p(0) 



(3) 



We find that the particle density decays fast in short time 
regime, followed by a slow decay in the long time regime 
as shown in FiglT] Indeed, numerically obtained values 
listed in Table [T] arc close to the one obtained from the 
formula c?s/4, and different from the ones obtained from 
the formula 
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FIG. 1: (Color on line) The particle density as a function of 
time on the critical branching trees with various degree expo- 
nents. Guidelines have slopes 0.40 (top) and 0.34 (bottom). 



Next, we study the reaction kinetics on deterministic 
fractal SF networks, introduced by Rozenfeld et al. [i^ . 
the so called (m, i))-flower and (u, t')-tree networks. These 
networks arc hierarchical networks, generated iteratively 
from a simple basic structure to higher level ones. Each 



TABLE I: Comparison of the exponent a numerically ob- 
tained, denoted as Onum, with ds/4 for various degree ex- 
ponent 7's for the critical branching tree. For comparison, we 
also list the mean field value obtained from the formula ([T]). 
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Q^num 


d,/4 


MP value 


2.5 


0.4 


0.38 


2.00 


2.7 


0.36 


0.35 


1.43 


3.5 


0.35 


0.33 


1.00 


4.0 


0.34 


0.33 


1.00 


4.5 


0.34 


0.33 


1.00 



link in the ?i-th generation is replaced by two parallel 
paths of u and v links long. Detailed rule can be found in 
Ref. (2^ . Depending on the rule, constructed networks 
are either the flower structure which contains loops or 
trees. These networks are fractal SF networks with the 
degree exponent, 7 = 1 -f ^^^^tj^, the fractal dimension, 
dt = ^"^"+") and the spectral dimension, = ^^"("+") 
for flowers, and ^ for trees. Numerical values of 
the exponent a are close to those from a = c?s/4 as can be 
seen in Tables II and III for the flower and tree structures, 
respectively. 
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FIG. 2: (Color on line) The particle density oi Aot B species 
as a function of time for the [u, ?;)-flower networks. Guidelines 
have slopes 0.43 (top) and 0.38 (bottom). 

To see if the segregation of yl-rich or i?-rich domains 
forms, we examine a quantity. 



QAB{t) = 



Nab 



Naa + Nbb ' 



(4) 



where NAB{t) is the number of {A,B) pairs located at 
the nearest neighbors averaged over different initial con- 
figurations. Naa and Nbb are similarly defined 23|. If 
Qab 0, then there is few pairs of different species at 
neighbor nodes, whereas if Qab 1, particles are mixed 
randomly. Since the particle density decreases in time. 
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TABLE II: Comparison of the exponent a numerically ob- 
tained, denoted as Cfnum, with ds/4 for various degree expo- 
nent 7's for the (u, ?;)-flower networks. For comparison, we 
also list the mean field value obtained from 



{u,v) 
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Q^num 


d./4 


MF value 


(2,2) 
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0.53 


0.5 


1.0 


(2,4) 


3.58 


0.45 


0.43 


1.0 


(3,3) 


3.58 


0.43 


0.41 


1.0 


(2,6) 
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0.43 


0.42 


1.0 


(4,4) 
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0.38 


0.38 


1.0 
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TABLE III: The same as Table II for the (u, v)-tree networks. 



(u,v) 
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f^num 


rfs/4 


MF value 


(2,2) 
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0.34 


0.33 


1.0 


(2,4) 


3.58 


0.37 


0.36 


1.0 


(3,3) 


3.58 


0.31 


0.31 


1.0 


(2,6) 
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0.38 


0.38 


1.0 


(4,4) 
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0.31 


0.30 


1.0 



FIG. 3: (Color on line) Plot of Qab, Nab and Naa versus 
time t for a (3,3)-flower network. The slopes of guidelines are 
-0.65, -0.78 and -1.43 from the top. 



TABLE IV: Comparison of the exponents for Naa and Nab 
between theoretical and numerical values for various em- 
bedded spaces, one dimensional regular lattice (1 dim), two 
dimensional square lattice (2 dim), critical branching trees 
(CBT) with 7 = 2.5 and 7 = 4.0, and (3,3)-flower hierarchi- 
cal network. 



their separation becomes large and two particles hardly 
locate at the nearest neighbors. We examine Nab and 
Naa independently as a function of time. Interestingly, 
they decrease with time in a power-law manner as shown 
in Fig. [3l which can be explained as follows: 

First, we examine Naa- The linear size £d of a do- 
main containing a species grows with time as t^/du, ^ 
A typical closest distance iAA between two particles of 
the same species scales as {l/pY^'^f. Assuming that 
p{t) one can obtain that £aa ^ ti/(2d») 

When ds < 2, the case of concern in this Letter, random 
walks are compact within the diffusion volume , and 

d f 

thus that is also valid within the volume (-aa- "^^^ prob- 
ability to find two such particles at the nearest neighbors 
is l/t'^iA- Thus Naa scales as {I / t/j^) p{t) . That is, 

NAA{t)^r''^'^. (5) 

Second, we examine NAsit)- When two particles of 
different species arrive at the nearest neighbors in the 
diffusion process, they can annihilate at the next step 
with a finite probability. Thus, we may set NAsit) oc 
dp/dt, and obtain that 

NAB{t)-t-''^/^-\ (6) 

Next, Qab is obtained as Nab /Naa- We compare the 
results obtained from simple arguments with numerical 
ones in Table Hvl 

To confirm that the segregation is caused by local hubs 
in the fractal structures, we destruct the local hubs by 



space 


N, 


lA 


Nab 






ds/2 


Num. 


(rfs/4) + 1 


Num. 


1 dim 


0.5 


0.49 


1.25 


1.27 


2 dim 


1.0 


0.99 


1.50 


1.59 


CBT (7 = 2.5) 


0.75 


0.59 


1.38 


1.28 


CBT (7 = 4.0) 


0.67 


0.63 


1.33 


1.28 


(3,3)-flower 


0.82 


0.78 


1.41 


1.43 



rewiring the links in the (3,3)-flower network while con- 
serving the degree distribution. Fig. 2] shows that the 
exponent a changes from a « 0.43 to the mean field 
value a = 1 as the number of rewired links increases. 
Moreover, Qab docs not decrease monotonically for the 
rewired networks as shown in Fig. [51 

While many complex networks in real world are frac- 
tals, the Internet at the autonomous system level is not a 
fractal. This may be caused from the geographical effect. 
Due to this non-fractality, the segregation does not occur 
in the Internet in the two-species annihilation, and thus 
the particle density decreases fast with exponent a w 1.8 
from recent Internet topology in the year 2004 as shown 
in Figini This property can be used beneficially when 
one designs a protocol for P2P network, virus-antivirus 
annihilation robot, etc [j^ . 

It is noteworthy that whereas the decaying behavior 
obeying the formula llj applies to the BA model when m, 
the number of incoming links at each time step, is larger 
than 1, it is not so for the BA tree network with m ~ 
1. This is because the tree structure has limited paths, 
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FIG. 4: (Color on line) The particle density versus time on 
rewired networks from a (3,3)-flower network. The slope in- 
creases as the fraction / of rewired links increases from / = 
to / = 0.5. For the / = and / = 0.5 cases, the slopes are 
close to a ~ 0.43 and 1, respectively. Since 7 « 3.58 > 3, 
a = 1 is the mean-field result. 
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FIG. 5: (Color on line) Plot of Qab as a function of t for 
the rewired networks used in Fig. U) Slope of the solid line is 
-0.65. 
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FIG. 6: (Color on line) The particle density as a function of 
time on the Internet in several difi'erent network profiles, the 
years 1997, 2000, and 2004. The slope of guideline is 1.8. 



which enhances segregation. Thus, the particle density 
decays slowly with exponent a « 0.5, even though the 
BA tree is not fractal. 

In summary, the segregation effect in the two-species 
annihilation reaction dynamics has to be taken into ac- 
count when the dynamics takes place on fractal, modular, 
or tree networks. In this case, the role of hubs is different 
from that of random SF networks and the particle den- 
sity decays slowly in a power-law manner with exponent 
less than 1, even though those networks are scale free. 
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